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ABSTRACT 

Context. The determination of cluster masses is a complex problem that would be aided by information 
about the cluster shape and orientation (along the line-of-sight). 

Aims. It is in this context, that we have developed a scheme for identifying the intrinsic morphology and 
inclination of a cluster, by looking for the signature of the true cluster characteristics in the inter-comparison 
of the different deprojected emissivity profiles (that all project to the same X-ray brightness distribution) and 
by using SZe data when available. 

Methods. We deproject the cluster X-ray surface brightness profile under the assumptions four different 
ellipsoidal geometry and inclination configurations that correspond to four extreme scenarios, by the non- 
parametric algorithm DOPING. The formalism is tested with model clusters and then is applied to a sample 
of 24 clusters. While the shape determination is possible by implementing the X-ray brightness alone, the 
estimation of the inclination with the line-of-sight is usually markedly improved upon by the usage of SZe 
data that is available for the considered sample. 

Results. We spot 8 prolate systems, 1 oblate and 15 of the clusters in our sample as triaxial. In fact, for 
systems identified as triaxial, we are able to discern how the three axis ratios compare with each other. This, 
when compounded by the information about the line-of-sight extent, allows us to constrain the inclination 
quite tightly. 

Key words. Methods: analytical; Galaxies: clusters: general 



1. Introduction 

The identification of the three dimensional cluster 
shape is of significant importance in the pursuit of 
quantities of cosmological interest. In particular, the 
misinterpretation of the true cluster geometry can af- 
fect the extraction of the Hubble constant from X- 
ray and SZe information (IZaroubi et al.l I2001I) and 
of the cluster mass. Erroneous cluster masses are of- 
ten the cause of less reliable understanding of cluster 
physics and poorer constraints from cluster cosmol- 
ogy. A mandatory exercise that should be undertaken 
towards the extraction of the correct cluster masses 
is the identification of the correct geometty and in- 
clination of the cluster. However, the lack of knowl- 
edge about these attributes are typically bypassed by 
resorting to the assumption of sphericity, even when 
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the projected ellipticity of the sy stem is indicated to 
be non-zero from obser vations dFabian et alJll98lt 
lYoshikawa & Sutolll999l;|Pizzolato et alj|2003l) . This 
is also the adopted geometry in frequently used pack- 
ages for determini ng dark matter distr ibutions in 
clusters, like JACO (Ma hdavi et~ai1l2007l) . 

On the contrary, halos realised in cosmological 
simulations have been found to be flattened and triax- 
ial in shape; the distribution of ellipsoid al shapes of 
simula t ed halos have been dis c ussed b y iFrenk et alJ 
dl988l) ; iDubinski & Carlberd (fl99l; IWarren etaU 
dl992|) : ICole&Laceyl dl996h : iBailin & Steinmetzl 
( 120050 , among others. Moreover, the formation of 
large scale structure has a be aring on the intrin - 
sic shape of cluster sized halos (Plionis et al. 2004); 
thus, the 3D cluster morphology could be a device 
that can be used to constrain cosmological models. 
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While the statistical analysis of th e 3D cluster 
morp h ology has been looked i nto (see iJing & Sutol 
l2002h iKasun & Evrardl 120051; lHavashi et alj 120071 
and references therein), the determination of the full 
3-D shapes and orientations of individual halos is a 
harder problem that has a ttracted relatively less at- 
tention (ISereno et al.ll2006l) . 

Recovery of the full 3-D morphology of a cluster, 
using measurements of accessible cluster character- 
istics (such as X-ray emission profiles or 2-D X-ray 
surface brightness distributions or SZe data), would 
require expertise in deprojecting such observed in- 
formation into the full 3-D distribution, under gen- 
eral triaxial geometries. Such can be achieved via 
parametric fits; in fact, parametric deprojection of the 
observed optical surface brightness maps of galaxies, 
under assumptions of sphericity and (more rarely) 
axisymmetry, ha s been studied before ( Palmer 1994; 
IB endmelhl 1 1 99 1 1) . However, the fundamental prob- 
lems with parametric fits are (i) the answer depends 
on the choice of the parametrisation and (ii) the 
goodness-of-fit quantifier (such as a x 1 measure) 
can appear spuriously inflated, particularly in the 
presence of non-homog eneous measurement noise 
dBissantz & Munkll200lb . 

Thus, non-parametric deprojection is a better op- 
tion. However, deprojection algorithms that promise 
improved three dimensional cluster mass distribu- 
tions, by taking the measured ellipticity and the true 
cluster morphology a nd orientation i nto ac count, are 
limited in availability. IZaroubi et alJd200ll) report on 
the application of a non-parametric deprojection al- 
gorithm that assumes axisymmetry, to a set of simu- 
lated clusters. However, as the authors state, the ap- 
plicability of this scheme to the current state of data 
sets appears difficult. 

A comprehensive approach that calls for the 
amalgamation of two or more mass indicators in- 
clude the implementation of X-ray observations, SZe 
information, lensing results and dynamical measure- 
ments. Such exercise s have been undertaken already 
(IZaroubi et alj Il998t iReblinskvl l2000t ISand et alj 
l2002f). An attempt has been made in the recent past 
dSereno et alj 120061) to decipher the three dimen- 
sional morphology and inclination distr ibution of a 
sample of clusters from Reese et al. (2002), with the 
aim of improving the cluster mass estimation (from 
X-ray measurements) more accurately. This work 
emplo ys the rudimentary SZe temperature decrement 
data dDe Filippis et al.l 120051) in the /^-modelling of 
the cluster. Thus, this work is susceptible to very 
large error bars that currently plague SZe data. This 
is of course topped by the errors that can be intro- 
duced by the choice of ellipsoids of revolution that 
are implemented to model their sample clusters. 

In this work, we present a novel, model- 
independent trick to determine the correct morphol- 
ogy and inclination of a cluster, without resorting 
to the assumption of axisymmetry. Such a deter- 



mination is made possible by analysing the multi- 
ple deprojected emissivity or X-ray luminosity den- 
sity distributions that are recovered under distinct 
assumptions about the cluster geometry and incli- 
nation and that will all project to the same ob- 
served X-ray brightness map. Thus, the deprojection 
in question needs to be performed by a scheme that 
is efficient in carrying out deprojection under gen- 
eral geometries, at any assumed inclination. This is 
possible with a new inverse deprojection algorithm 
DOPING (Deprojection of Observed Photometry us- 
ing an INverse Gambit ) that has been reported in 
IChakrabarty & Ferraesd (120081) . 

The inter-comparison of the amplitudes and 
shapes of the different deprojected emissivity distri- 
butions tell us the correct shape of the cluster (includ- 
ing triaxial ones) while the inclination constraints are 
improved upon by the implementation of the avail- 
able cluster elongation information from the SZe 
data, as reported in the literature. In case of triaxial 
systems, the SZe data is used to determine the two 
intrinsic axial ratios. 

The paper is organised as follows. Following the 
first introductory section, we discuss the deprojec- 
tion algorithm in brief. Section 3 is devoted to the 
methodology that we use to extract the cluster shapes 
from the deprojected profiles. The testing of the ad- 
vanced scheme is dealt within the following section. 
We then proceed to apply this method to our cluster 
sample (Section 5), which precedes the results sec- 
tion. The paper is rounded off by recounting some 
relevant aspects of the presented work. 



2. DOPING 

DOPING stands for Deprojection of Observed 
Photometry using an INverse Gambit. However, the 
name does not do full justice to the capacity of the al- 
gorithm; it is a deprojection algorithm that could just 
as well be applied to deproject observed X-ray 2-D 
surface brightness maps (SB) of clusters. The depro- 
jection can be performed in general triaxial geome- 
tries and is able to incorporate radial variations in 
shape. The deprojected luminosity density distribu- 
tion is sought by penalised likelihood approach that 
uses an MCMC optimiser. Although DOPING is in 
general able to use the 2-D S B information and pro- 
vide the full three dimensional p, in the case of an 
oblate/prolate system with a uniform projected ec- 
centricity, the effective dimensionality of the prob- 
lem reduces by one, since the projected and intrin- 
sic eccentricities can then be related through analyti- 
cal relations. This idea is exploited in the work. The 
S B is treated as a function of distance along a photo- 
metric axis (the x-axis). Similarly, the emissivity or 
X-ray luminosity density distributions are presented 
along the x-axis. 
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3. Method 

As explained in Section[2] DOPING is able to incor- 
porate information about the LOS extent of a clus- 
ter from SZe measurements, in the process of depro- 
jecting under the assumption of triaxiality. However, 
the poor quality of the currently available SZe data 
will cause the determination of shape ambiguous. 
Thus, for the current purpose, we will be carrying 
out the deprojection of the observed X-ray SB un- 
der assumptions of prolateness or oblateness for the 
intrinsic shape of the observed cluster; as for its in- 
clination, we choose to deproject at the two ends of 
the range of inclinations allowed for an oblate geom- 
etry (discussed below). The other, even more impor- 
tant reason for choosing to deproject at these geome- 
try +inclination configurations is the relative simplic- 
ity of deprojection under oblateness/prolateness than 
triaxiality. 

The different X-ray luminosity density distribu- 
tions or the emissivity distributions (p) that will be 
recovered from such deprojections, will be analysed 
to tell us if the system under consideration is prolate, 
oblate or even triaxial. 

Now, let us examine our choice of inclinations 
for deprojection, in greater detail. We will carry out 
the deprojections under the assumptions of i-90° and 
i=imin, where z' mI „ is the smallest allowed inclination 
for the measured projected eccentricity for the as- 
sumption of oblateness (/„„■„ = sin -1 e p , where e p is 
the projected eccentricity that is related to the pro- 
jected axial ratio q p as e 2 p — 1 - l/q'i). As a result 
of such deprojection of an S B, we will be able to say 
if the cluster is inclined to the line-of-sight (LOS) by 
an angle that is less than /„„■„ or is of a value interme- 
diate to 90° and /„„■„. 

Thus, there are 4 distinct deprojection configura- 
tions (D) that we will consider for any given SB: 

1. assumed geometry is prolateness and assumed 
i=90°: D p 90 , 

2. prolateness and z'=z mI „: D'' nin , 

3. oblateness and /=90°: D^ Q and 

4. oblateness and i = /„,„: D" . . 

mm mln 

On the basis of these deprojections, we will identify 
the shape and inclination of a given cluster. 

3.0.1. Effect of D 

It is to be emphasised that our determination of the 
true 3-D shape of a cluster is not affected by the exact 
geometry +inclination configuration that we choose 
to deproject its X-ray S B under; rather, it is the inter- 
comparison of the various profiles of p that are recov- 
ered from deprojection under the different geome- 
try+inclination configurations that we choose. Thus, 
if we choose to deproject an observed X-ray SB, as- 
suming the cluster to have a geometry +inclination 
different from the ones that we choose, then the inter- 



comparison of the corresponding profiles of p will 
be different from what we infer here. It is a question 
of standardising such an inter-comparison, in order 
to be able to conclude the true cluster characteris- 
tics. All we do in this paper is advance the same, us- 
ing the trends in the inter-comparison that we notice, 
when the cluster at hand is deprojected under our 
choices of geometry +inclination, i.e. D. The modus- 
operandi that connects such an inter-comparison to 
the true cluster characteristics, is discussed below 
(SectionUSll. 

3.1. Coordinate System 

The two coordinate systems that suggest themselves 
readily to any projection related investigation are the 
body coordinates of the observed system (considered 
regular Cartesian, marked here in upper case letters 
X — Y — Z) and the observer's coordinate system 
(marked in lower case x—y — z). The POS is spanned 
by the x—y plane while z is along the LOS. Also, we 
consider one of the principal axes of the cluster (X) 
to be coincident with one of the photometric axes x 
(the photometric major axis for oblate models and 
the photometric minor axis for the prolate case). 

3.2. Underlying Principle 

The determination of a cluster's intrinsic shape 
and inclination is based on the fact that the inter- 
comparison of p recovered by deprojecting under the 
assorted geometry and inclination configurations D, 
bears the signature of the true characteristics of the 
cluster. 

The comparison of (the shape and amplitude of) 
two distinct profiles of p that will be recovered from 
deprojecting a given X-ray S B, under two distinct de- 
projection scenarios (2 different Ds), will depend on 
the inter-comparison of the 

- extents along the LOS, of the 2 systems that are 
ascribed to the 2 different assumptions of geom- 
etry (corresponding to the 2 Ds) 

- ellipsoidal radii to the same generic point (x, y, z), 
in the 2 systems described by the 2 Ds, assuming 
the gas density to be stratified on concentric el- 
lipsoidal shells. 

In this section, these ideas are expounded upon. 

3.2.1 . Effect of Varying the LOS Extent 

The luminosity density that is required to be pro- 
jected through a given length, in order to produce an 
observed S B, must be higher when the extent along 
the axis of projection is less in a deprojection model 
than in the case of the "true system", i.e. the observed 
cluster. In other words, LOS d < LOSj => Pd > Pt- 
Here, the subscript "D" refers to quantities relevant 
to the deprojection model while the subscript "T" 
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refers to quantities relevant to the true system. (We 
will follow this notation elsewhere in the paper too). 

3.2.2. Effect of varying the 3-D Density Structure 

We can develop a qualitative understanding of the 
shapes of the recovered density profiles by assuming 
that the density is given as a function of the ellip- 
soidal radius The form of £ is determined by the 
geometry and inclination i that describe a system, as 
well as its intrinsic eccentricity e. Thus, for an oblate 
system, inclined at angle i, for any (x, y, z), £ is: 

, , , sin 2 ; 9 9 cos 2 / 

x L + y z (cos L i + T ) +z z (sin z / + .-) (1) 

1 - e L 1 - e L 

■ T ^ 

-yzsin2i- j = £ . 

1 - e z 

Thus, we assume p = p(£). We also assume that the 
physically motivated idea of monotonicity holds true 
in the density distribution: p(£,-) < p(£j), V£; > tjj. At 
this stage, this exercise is held independent of effects 
of LOS extent on the density. 

Now, let us consider two systems with different 
e, i and geometries. 

- Let the systems be such that in System I, the el- 
lipsoidal radius to a general point x, y, z is higher 
than in System II, i.e. to x,y,z, > (say) 

- Then pi(x,y,z) < pil(x,y,z), since we assume 
density to monotonically fall with 

- Since the projection of density (observed S B) is 
given from measurement, to compensate for this 
trend, in System I, a given value of density would 
show up at a higher than in System II. 

This translates to the given density showing up at a 
higher x, in System I than System II, when the den- 
sity profile along the x-direction is sought. 

Now let us replace System II with the true system 
and let System I be one of the 4 deprojection models 
that we consider. Thus, in this case, go > £r- This 
comparison can imply stronger effects on the density 
at smaller values of x or larger values of x, depend- 
ing on the relative geometries of the two systems in 
question. While such situations are explored below, 
now let us monitor the consequences of £o exceed- 
ing first near the core of the cluster and then at 
the outer parts of the cluster. 

When £d > %t around the core, it has the fol- 
lowing effect on the shape of the recovered density 
profile. If the density profile of the true system is as 
shown in solid lines on the left in Figure |2j then the 
recovered density profile will be flatter at the centre, 
as represented by the broken lines in this figure. If 
the smallest x at which density information is given 
for the true system is xq, then that for the deprojec- 
tion model will be x\, where x\ > xq, for a given 
choice of the radial binning. In lieu of any informa- 
tion available for the radial bins inner to X\, the code 




- 11 - 

Fig. 1. The ellipsoidal radius to a point P is marked as £1 
and £2, for the two different ellipsoidal configurations, the 
projection on the plane of the paper of which are shown 
in thick solid and thick broken lines. The system with the 
bigger long axis and the smaller short axis corresponds to 
a higher ellipsoidal radius. Along with the variations in 
depth, differences in ellipsoidal radius is the other factor 
that is affected by the choice of the deprojection model. 

will spread light out uniformly over these bins, using 
the same density that is recovered at x\. Thus, the 
deprojected density profile will appear to manifest a 
larger core than the true density profile. 

However, when the inequality: £d > %j is 
stronger at larger than at smaller x, the deprojected 
density profile will end up being flatter on the out- 
side than the true system. 

In contrary to this case, if we have a deprojec- 
tion model that generally corresponds to a lower ^ 
at any x,y on the POS, than in the true system, the 
code will attribute a given density to a lower value of 
x in the density profile recovered with the deprojec- 
tion model than in the true system. The deprojected 
density profile will then appear to manifest a similar 
shape to the true density profile, but will have a com- 
paratively lower amplitude. This case is shown in the 
right panel of Figure|2] 

In general, we can base our understanding of 
the density structure in a cluster, by relying on the 
following ideas that are illustrated in Figure [2] and 
Figure [T] 

- In two prolate systems with the same long axis, 
the cluster with the relatively smaller short axes 
will imply a bigger £ at a general point, as long 
as the point in question lies close to the long axis 
of the system. 

- Similarly, in two prolate systems with the same 
short axes, the bigger is the long axis, bigger is £ 
to a general point, as long as this comparison is 
sought within the extent of the smaller of the two 
long axes in question. 
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- In two oblate systems with the same long axes, 
the smaller the short axis, the bigger is ^ to a gen- 
eral point x, y, z, (the point is considered a part of 
both systems). 

- Similarly, when we consider two oblate systems 
with the same short axes, the cluster with the big- 
ger long axes corresponds to a bigger £ . 

- In general, it is true that a prolate system with the 
same axial ratio as an oblate system, will imply a 
larger £ to the point x, y, z than the oblate system. 

- It is also generally true that deviations of an ob- 
served cluster from a prolate deprojection model 
will affect the outer parts of the density profiles 
while differences between an observed system 
and an oblate model implies differences in the re- 
covered density profile that will manifest them- 
selves at lower values of x. 

3.2.3. True system: T? 

Let us elucidate the effects discussed in Section [3.2.1l 
and Section [3.2.21 in reference to the true system T P , 
i.e. a cluster that is prolate in shape and is inclined 
at a angle i\i m j„ < i < 90°. The projected axial ratio 
is q and the extent along the photometric semi-major 
axis is a. 

A generic T? system is shown on the left in 
Figure 21 the X-ray SB observed from this true sys- 
tem is then deprojected under the 4 deprojection sce- 
narios that are shown in progression in Figure [4] In 
the context of this figure, the truly prolate system is 
projected along an inclination such that the projected 
axes on the POS are b (< a) and a/q (by construction 
of the coordinate systems). If we now assume that 
this projection on the POS is due to a system that 
is prolate and inclined at 90° to the LOS, then that 
system will be ascribed intrinsic axes of a/q,a/q&b. 
If however, the system is assumed to be prolate and 
inclined at an angle /„„■„ (< /), the long axis of the 
system will have to exceed that of the true system, in 
order for this configuration to project to the observa- 
tions. The state of the system, if assumed oblate and 
inclined at either 90° or i mi „, is shown in the second 
from right and right-most positions in Figure |4] 

Thus, we see that given a 2-D image, it is possi- 
ble to guess the semi-axes lengths of the 4 systems 
that can be described by the 4 different D that we use 
herein. Now the question arises as to how the run of 
p with x is dictated by the lengths of the semi-axes in 
the corresponding configuration; more to the point, 
we try to gauge how p recovered for a given D will 
compare to the that for another choice of D, as well 
as the true system. The key to such trends lies in the 
mutual weighing of two configurations for their LOS 
extents and the ellipsoidal radii to a generic point. 

For instance, given that the long axis of D^ Q is 
shorter than that of Tf, we realise that the ellipsoidal 
radius to a point is higher in T p than in Dg Q (2 nd 
itemised point in the last list in Section [3.2.2t . Thus, 



the density profile recovered by deprojecting under 
Dq will have a similar shape as the true p, but will 
be lower in amplitude. However, the details of the 3 
semi-axes lengths of the D'^ suggests that the LOS 
extent is lower in D P Q than in T? , Thus, the lessen- 
ing of the amplitude caused by differences in £ is 
compensated by differences in LOS extent. Similarly, 
we can track how the p recovered for the other 3 
D will tally with the true density profile, as well 
as how these recovered profiles compare with each 
other. The details of the other cases are listed below. 
We then register such an intercomparison as possi- 
ble, when the true system is T p with i e (i m ;„,90°). 
Thus, whenever the S B of an observed system, when 
deprojected under the 4 used Ds, correspond to p pro- 
files that all have the same shape and the amplitudes 
of which tally in the way described below, we turn 
our experience around and declare the system to be 
prolate in nature with inclination < /„„■„. 

Had the inclination of the true system at hand 
(Tf) been such that i < z' m! „, then even for D p . , the 
long-axis would have fallen short of the long axis of 
the true system, as for the D p case. In this case, the 
recovered profile of p would have compared with the 
true density profile similar to the p recovered for Dg Q . 
Thus, the intercomparison of the density profiles re- 
covered in this case is registered as possible when the 
true system is Tf, with i e (0, /,„,„). 

We qualify the word "possible" in the two in- 
stances above, simply because we are finding a pos- 
sible answer to the question of: what is the correct 
shape and inclination of the cluster that projects to 
the observations. In other words, we do not claim 
that our answer is unique. It maybe possible in 
some circumstances, to conjure up a different geome- 
try+inclination specification, that will correspond to 
the same projection. We discuss this in Section|7] 

Had we used other deprojection configurations, 
the intrinsic axial ratios in such configurations would 
have been different from what they are for the 4 D 
that we use. The intercomparison of the profiles for p 
that emanate from such conjured configurations (the 
D of our choice), is rendered easier to interpret given 
that the D of our choice correspond to the extrema of 
the geometry and inclination scales. 

Some details: 

- intrinsic axial ratio of this true system is higher 
than the deprojection models inclined at 90°, 

- with the intrinsic major-axis exceeding that of the 
Dg case but the intrinsic minor-axis exceeding 
that of the D° min model; 

- the intrinsic major axis of the true system falls 
short of the same in the D p . but the intrinsic mi- 

min 

nor axis falls short of that in the D^ Q model. 

For the 4 different deprojection models, the char- 
acteristics are enumerated as follows: 
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Fig. 2. If the ellipsoidal radius to most points on the z = plane in the deprojection model is less than that in the true 
system, then density of the deprojection model will in general be higher than in the true system. In this case, the density 
profile recovered along the x-axis will be characterised by the same value of density appearing at smaller values of x. 
This results in a deprojected density profile that is as shown in broken lines on the right panel. The true density profile 
of the observed cluster is shown in the same panel, in solid lines. The deprojected profile is shown schematically as 
simply the amplitude-lessened version of the true profile and therefore is shown to extend to negative x-values; in reality, 
deprojected density values are recovered only for positive x. On the other hand, if g D exceeds g T on the average, then 
deprojected density will fall short of the true density, which will imply a larger core size in the profile recovered in this 
case (shown in broken lines on the left). 



- LOS d < LOS t => amplitude of po at the 
centre exceeds that of p T . 

- Ellipsoidal radius to a general point x, y, z is 

£z> < £,j =>Pd > Pt- 

- This results in deprojected density profile 
(p x D ) showing a similar shape as p x T , but is 
lower in amplitude. 

y. : 

mm 

- LOS D > LOS T ==> amplitude of p D falls 
short of pr- 

- h > fr. 

- The density comparison elicited by this in- 
equality implies that on the outside, p x D is 
lower in amplitude and is also flatter than p x T . 
This flatness will be less than that for a simi- 
larly eccentric, true oblate system, for this de- 
projection scenario (see next sub-section). 

- LOS d — a, which is in excess of LOS j => 
amplitude of p D falls short of that of p T . 

- £d < fr- 

- This implies that p x D shows up with a shape 
that is similar to the true system but is lower 
in amplitude; in fact the recovered profile is 
similar to that recovered from Lf min . 

y. : 

mm 

- LOS d < LOS t ==> amplitude of po > pj. 

- Given that this case corresponds to a smaller 
minor axis than the true system, and that the 
true system is described by a comparatively 



larger intrinsic major axis, %o « for lo- 
cations that overlap between the true system 
and deprojection model. 
- p A (Dg ) has similar shape as the true system 
and is highest in amplitude out of the 4 de- 
projection models. 



3.2.4. True system T° 

Following similar logic as used in the last section, we 
conclude that when the observed cluster is oblate, in- 
clined at an angle i\i min < i < 90°, the density profiles 
recovered from the 4 deprojection models are as fol- 
lows. 



- Dg Q : recovered density profile will have similar 
shape as the true system but will be higher in am- 
plitude from LOS considerations. 

- Lf min : recovered density profile will manifest a 
larger core than the true profile and higher in am- 
plitude from LOS considerations. 

- Dg Q : similar shape to true profile but amplitude 
will be lower from considerations of ellipsoidal 
radius but can be either lower or higher than the 
true profile, from LOS considerations. 

- D° min \ recovered profile will have the same shape 
and higher amplitude. 
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Fig. 3. The three kinds of triaxial systems encountered in this work - Type I and Type II (left, b > a > c) and Type III 
(right, b < a < c). The high inclination and low bja ratio in Type II systems ensures that the photometric major axis is 
along the j-direction, unlike in Type I systems. 



3.2.5. Triaxial Systems 

Here are some types of triaxial systems that we can 
identify via the inter-comparison of the deprojected 
luminosity density profiles in our sample clusters. 
This is discussed in terms of the extent along the 
principal axis that is identified as the LOS for an in- 
clination of 90° - referred to as the F-principal axis 
in the following discussion. 

- A triaxial configuration that deviates from a per- 
fectly oblate shape in that the F-principal axis is 
the longest of the three while the principal axis 
that corresponds to a photometric axis being the 
second longest - this is referred to as "triaxial- 
Type I" (left-most panel in Figure|3). 

- It is also possible that the true system is triaxial 
in the same way as Type I systems but the extent 
along the X-axis is only slightly smaller than that 
along Y, so that for inclinations close to 90°, the 
photometric semi-axes in the ^-direction exceeds 
that along the y-direction, in contrast to the POS 
picture for Type I systems. Such a triaxial config- 
uration is referred to as "triaxial-Type II" (middle 
panel in Figure [3}. 

- If the sole deviation from pure prolateness that 
is responsible for triaxiality is such that the Y- 
principal axis is the shortest of the three, while 
the principal axis on the POS has intermediate 
extent, then the ensuing triaxial configuration is 
referred to as "triaxial-type III" (right panel of 
Figure 0). 

- Since the relation between the projected and in- 
trinsic axial ratios in the triaxial case is differ- 
ent from that of the oblate and prolate cases, it is 
not possible to extract accurate constraints on the 



system inclination from the deprojection work 
alone. However, in lieu of any other information, 
we can make an approximate guess about what 
the inclination range for a certain triaxial cluster 
is. This is done by suggesting if the inclination is 
intermediate to 90° and /„„■„ or less than /„„,,. The 
correct estimate of the allowed inclination range 
for a triaxial cluster will be made on the basis of 
the pair of relations (that are discussed below), 
between the projected and intrinsic axial ratios 
of the triaxial system. 



3.2.6. Other Important Predictions 

We can pursue logic similar to what was used above, 
to arrive at the following, critically important conclu- 
sions. 



- When the true system is TL, all profiles will be 
similar in shape, with the relative amplitudes of 
the profiles recovered under the 4 different depro- 
jection scenarios, as for a truly prolate system, 
inclined at a general angle i. 

- When the true system is TS Q , all profiles will be 
similar, with the amplitudes as in the case of in- 
clination at a general angle i. 

- For a truly prolate cluster, the deprojected pro- 
files will be further apart, the smaller the angle 
of inclination is. On the other hand, for a truly 
oblate cluster, the amplitude differences between 
the recovered density profiles will increase with 
decreasing inclination. 



True - Prolate, i 

A 

a/q 

a/q 
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Fig. 4. Figure showing interpretation of the observed prediction of a cluster that is truly prolate and inclined at an angle 
i e [0, 90°] (shown on the top left-most corner), under 4 deprojection scenarios: Dg Q (second from left), D^. (third from 
left), D° go (second from right) and Djj. (right). The cluster is inclined such that its intrinsic axial ratio q results in the POS 
projected axial ratio of b : a/q. The intrinsic semi-axes lengths that need to be ascribed to the 4 different D configurations 
are: ^ Q < f r , ^' mjn > f r , £JJ < <f T and %° njn w £ T . The comparison in terms of the LOS extent is such that only for D p min , it is 
in excess of that in the true system; for the other 3 D, it falls short of that in the true system. In light of these comparisons, 
(according to the last bulletted list in Section 3.2.2), the recovered deprojected density profile along x, is shown in the 
lower panels (in solid lines), for the corresponding D, along with the true density profile, plotted in broken lines. 



4. Tests 

Here, we report the tests that were carried out to con- 
firm the predictions that were discussed in the last 
section. These tests were performed with 3 clusters 
that are either perfectly prolate or oblate, inclined at 
an angle intermediate to 90° and /„„■„, namely 60°. 
In order to demonstrate the independence of the re- 
sults from the core radius of the system, 2 of the 3 
test models were assigned a core radius of approxi- 
mately 0.01 arcmin and the other has a core radius 
of about 0.2 arcmin; all the test clusters are assigned 
a projected axial ratio of q p -l.2. Surface brightness 
data for each of these test models were estimated by 
performing a LOS integration of the X-ray luminos- 
ity density distributions that were chosen to admit an 
analytically integrable form. 

Each of the test brightness profiles was depro- 
jected under the 4 extreme deprojection scenarios: 
Do,,, D p . , D ( L and D° . The density profile that was 

90' mm' 90 mm J r 

recovered by deprojecting the brightness data of any 
test model, under any deprojection scenario, was re- 
projected back on the POS and this POS distribution 



was then compared to the input brightness data, to 
ensure the acceptability of the test run. 

The X-ray luminosity density profiles that were 
recovered with the different deprojection scenarios, 
for the 3 different test models, are displayed in 
Figure The similarity of the inter-comparison of 
the 4 different deprojected density profiles, in the 
cases of the two prolate clusters with the varying 
core radii, vindicates the core size independence of 
the shape determination. However, is this an accept- 
able conclusion, in terms of the expanse and details 
of the test systems we considered? 



4. 1 . Generality of Tests 

In this section we examine the fundamental ques- 
tion of how general our test systems are. To begin 
with, a statistical confirmation of our method is a far 
cry from the limited tests that we have presented be- 
low. It would be judicious to scan a range of ellip- 
ticities and core radii corresponding to the distribu- 
tion of these quantities, as constrained by observa- 
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tions and simulations. Furthermore, the redshift de- 
pendence of these distributions has to be folded in, 
to test for our shape determination technique (as in- 
dicated by the ellipticity-redshift relation). In fact, 
work is underway to calibrate our technique against 
a sample of simulated gas+dark matter halos from 
the Millennium Gas Simulations (Chakrabarty et. al, 
2008). 

While such a rigorous statistical test is on, we 
have presented here tests done with three clusters 
that have characteristics similar to what have been 
noted with the sample of clusters that we apply our 
method to (following section) - as apparent from 
Table 1, in our cluster sample, there are 7 out of 25 
clusters, that have a projected axial ratio lying be- 
tween 1.2 and 1.25, approximately (including A 399 
which is at a redshift of 0.072). Most of the other sys- 
tems manifest even more ec centric 2-D images. In 
fact, iDe Filippis et al.l {2005) state that the redshift- 
ellipticity relation is not manifest in this sample. 
Thus, our test configurations as motivated by mea- 
surement. 



5. Application 

We a pply our formalism to a sample of 25 clusters 
from lReese et alj d2002f) . This sample consists of 18 
X-ray selected clusters with z > 0.14, 6 > -15° and 
Lx(0.1-2.4keV) > 5xl0 44 /z 50 2 erg s" 1 andfor which 
high S/N detections of SZE, high-S/N X-ray imag- 
ing and electron temper atures were availab le. To 
these, 7 clusters from the Mas on et"aT1 d2001l) sam- 
ple were added, con taining clusters from X-ray flux- 
limited catalogue o fEbel ing et al.l(ll996h . Details on 
the completeness of the la tter subsample are given 
by iMason & Mversl d2000h . Actually, we choose to 
work with 24 of these 25 clusters since the measured 
surface brightness profile of one of the sample clus- 
ters (A 520) appears to be plagued by very large ob- 
servational error bars in the last 3 arc minutes (typi- 
cally, in excess of 50%). 

The motivation behind choosing this sample is 
that the three-dimensional intrinsic morphology of 
the clusters in this sample has already been re- 
ported, as estimated by an independent method 
that combi nes X-ray and Su nyaev-Zerdovich ob- 
servations (Sere no et al.l 12006). Thus, we will be in 
a position to com pare our results with that from 
ISereno et alJd2006l) . 

5.1. Data Analysis 

In our analysis we have used archival Chandra or 
XMM-Newton data for all clusters in the sample. The 
Chandra Interactive Analysis of Observation soft- 
ware (CIAO 3.3.0.1) and the XMM-Newton Science 
Analysis Software (SAS 7.0.0) were used through- 
out. 



We first modelled the projected two dimensional 
cluster emission, in the 0.3 - 7.0 keV energy band. 
Pixel values of all detected point sources were re- 
placed with values interpolated from the surrounding 
background regions. Using the SHERPA software, 
we fitted the cluster surface brightness to elliptical 
2D /?-models. The resulting best fit parameters for 
all clusters in the sample, relevant for this paper, are 
listed in Table Q] 

We then extracted one dimensional surface 
brightness profiles within elliptical concentric an- 
nuli, with projected ellipticity, position angle and 
centre matching those estimated in the our 2D analy- 
sis (as listed in TableQ~|i. The distance x from the clus- 
ter centre is measured along the major axis of the el- 
lipses. The profiles were extracted in the 0.3 -7.0 keV 
energy range. In the same energy band, the back- 
ground for each cluster was estimated, within a pe- 
ripheral region of the field of view. Corresponding 
radial profiles, extracted within the same elliptical re- 
gions from the exposure maps, provided the required 
average values of the effective exposure time within 
each annulus. Background subtracted and exposure 
corrected surface brightness profiles were then ex- 
tracted for each cluster. 

6. Results 

The intrinsic shapes and inclination classes that we 
infer for our chosen sample of clusters are enumer- 
ated in Table 2. We also include the sha pes and incli- 
nation s of the clusters, as predicted by ISereno et al.l 
(2006), in adjacent columns. 

Thus, we find that in our sample of 24 clusters, 
15 are triaxial, 1 is oblate and the rest (8) are pro- 
late. This shape distributi on is significantly d ifferent 
from that deciphered by ISereno et al.l d2006l) . as ev- 
ident in Table 2; we discuss the differences in detail 
in Section |7TT1 

We also looked at the central flux distribution 
for the clusters in our sample and sought correla- 
tions with the recovered shapes; in particular, we at- 
tempted to establish is any correlation exists between 
the clusters and the central flux measurements. This 
was found to be absent. 

6.1. Incorporation of SZe Data to Prolate & 
Oblate Clusters 

The impending improvement in the quality of 
the SZe data will boost this technique greatly. 
Knowledge of the extent along the LOS will per- 
mit us to perform the deprojection in a fully tri- 
axial geometry - such a deprojection is possible 
with DOPING. Thus, while in the current version 
of DOPING we compare the plane-of-the-sky pro- 
jection of the recovered density distribution to the 
observed brightness map of the system, on improve- 
ment of the SZe data, the LOS extent of the density 
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Fig. 5. X-ray luminosity density profiles of different test clusters, recovered by deprojecting the toy X-ray surface bright- 
ness maps of these systems. The test systems include prolate (right and middle) and oblate (left) models with core size of 
about 0.003 arcmin (left and right) and 0.2 arcmin (middle); the projected axial ratio is 1.2 for all three test systems. The 
density profiles recovered under the D^ Q scenario is in magenta, the D p mjn case in blue, D° 90 in black and D° ljn in green. 



structure can also be pinned down, for a given choice 
of the inclination. We would then be in a position to 
better constrain both the intrinsic axial ratios for the 
triaxial case and thereby improve upon the estimate 
of the inclination of the cluster to the LOS. 

As mentioned above, one motivation for choos- 
ing this sample is the availability of the SZe data for 
all the clusters in the sa mple. This was reported by 
iDe Filippis et al.l d2005l) . It is not surprising that the 
ratio of the LOS extent to the photometric major axis 
(iws) °f a sample cluster is characterised by large 
error bars. 

We implement the reported q LO s first for all the 
clusters that are recovered as prolate in our sample, 
and then for the oblate o nes. To do th i s, we use the 
following relations (from lSereno et aL 

Jq 2 p - cos 2 i 

qws = — prolate (2) 

q p sin i 

yjl — q 2 p cos 2 i 

- ; oblate 

sin; 

The reported q p and qios values are used in this 
equation, for a sample cluster that has been identified 
as prolate. The values of q L0S reported at the two ex- 
trema of the error bands, for any cluster, results in a 
range for the inclination. These inferred range on the 
inclination is superimposed on the range of inclina- 
tions that were judged admissible from deprojection 
considerations alone. The range corresponding to the 
overlap of these two sets of constraints is then nar- 
rower, as represented in the last column of Table 3. 
A similar exercise was undertaken for the sole oblate 
cluster in our sample (according to our shape evalua- 



tion); the improved inclination estimate for this sys- 
tem is included as the last entry in Table 3. 

As can be appreciated from this table, the in- 
clusion of SZe data makes a significant effect in 
constraining the cluster inclination. However, the al- 
lowed inclination ranges from the X-ray data and the 
SZe data did not overlap for the case of Abell 611; 
thus, we are unable to produce an inclination for this 
cluster. 

6.2. Incorporation of SZe Data to Triaxial Clusters 

The direct exploitation of the SZe data is possible 
only for the cases of the purely oblate or purely pro- 
late clusters. For the general triaxial clusters, the SZe 
information is used to process the two intrinsic ax- 
ial ratios, given that we know how the lengths of the 
three semi-axes compare, given the recovered triaxi- 
ality type. Additionally, the cluster inclination can be 
pinned down from the recovered inequalities that re- 
late the intrinsic axial ratios to unity. Such constraints 
too are derived from the identification of the type of 
triaxiality in question. The details of the methodol- 
ogy used to constrain axial ratios and inclinations is 
delineated in Appendix A. 

As mentioned above, the inclinations guessed 
from the X-ray deprojection exercise alone are loose 
approximations only and only the values of i derived 
from the incorporation of the SZe data should be con- 
sidered correct. 

The axial ratios in the considered triaxial cases 
are tabulated, along with the improved inclination 
constraints, in Table 4. It is to be noted that these 
constraints a re estimated by ignoring the errors that 
are reported dDe Filippis et al.l[2005h on the elonga- 
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Fig. 6. X-ray luminosity density profiles of 6 of the 24 clusters in our sample. The colour scheme used for the profiles 
recovered under the 4 different deprojection scenarios is the same as in Figure|5] 
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Table 1. Cluster Sample 







z 




RA. 




Decl. 






9 






9 

deg 




Satellite 


MS 1137.5+6625 





784 


11 


40 


22 


3 


+66 


08 


15 


3 


1.113 


± 





014 


116.1 i 


: 1.0 


1 


MS 0451.6-0305 





550 


04 


54 


1 1 


4 


-03 


00 


51 


3 


1.307 


± 





015 


95.9 ± 


1.1 


1 


CI 0016+1609 





546 


00 


18 


33 


5 


+ 16 


26 


12 


9 


1.205 


+ 





013 


310.8 4 


: 1.7 




RXJ1347.5-1145 





451 


13 


47 


30 


7 


-11 


45 


09 


1 


1.453 


+ 





019 


158.5 4 


: 1.0 


1 


A 370 





374 


02 


39 


53 


3 


-01 


34 


39 





1.564 


± 





018 


6.2 + 


0.7 


1 


MS 1358.4+6245 





327 


13 


59 


50 


7 


+62 


31 


04 


1 


1.325 


+ 





019 


156.6 4 


: 1.4 


1 


A 1995 





322 


14 


52 


57 


9 


+58 


02 


55 


8 


1.242 


+ 





010 


57.8 ± 


1.0 


1 


A611 





288 


08 


00 


56 


8 


+36 


03 


23 


5 


1.14 


+ 





05 


34 ± 


9 


1 


A 697 





282 


08 


42 


57 


6 


+36 


21 


56 


8 


1.334 


± 





016 


163.8 4 


: 1.2 


1 


A 1835 





252 


14 


01 


02 





+02 


52 


42 


9 


1.225 


± 





012 


173.0 4 


: 1.4 


1 


A 2261 





224 


17 


22 


27 


1 


+32 


07 


57 


4 


1.022 


+ 





017 


90.0 ± 


1.7 


1 


A 773 





216 


09 


17 


53 


1 


+51 


43 


37 


9 


1.237 


+ 





022 


90.0 ± 


2.5 


1 


A 2163 





202 


16 


15 


46 


6 


-06 


08 


44 


9 


1.206 


± 





004 


90.0 + 


0.6 


1 


A 1689 





183 


13 


11 


29 


6 


-01 


20 


28 





1.141 


± 





012 


17.6 + 


2.2 


1 


A 665 





182 


08 


30 


57 


1 


+65 


51 


01 


8 


1.238 


+ 





012 


146.3 J 


: 1.2 


1 


A 2218 





171 


16 


35 


51 


9 


+66 


12 


34 


6 


1.162 


+ 





009 


96.5 ± 


1.5 


1 


A 1413 





142 


11 


55 


17 


9 


+23 


24 


16 


2 


1.473 


± 





019 


177.8 4 


:0.9 




A 2142 





091 


15 


58 


20 


1 


+27 


14 


03 


5 


1.540 


± 





007 


127.7 4 


:0.3 




A 478 





088 


04 


13 


25 


3 


+ 10 


27 


53 


5 


1.477 


+ 





006 


43.5 ± 


0.3 




A 1651 





084 


12 


59 


21 


9 


-04 


11 


44 


6 


1.184 


+ 





013 


87.6 ± 


1.9 




A 401 





074 


02 


58 


57 


1 


+ 13 


34 


37 


8 


1.303 


± 





008 


34.9 ± 


0.6 




A 399 





072 


02 


57 


52 





+ 13 


02 


38 


7 


1.207 


± 





009 


22.5 ± 


1.2 


2 


A 2256 





058 


17 


04 


00 


4 


+78 


38 


37 


1 


1.327 


+ 





008 


28.7 ± 


0.5 


2 


A 1656 





023 


12 


59 


44 


1 


+27 


56 


43 





1.141 


+ 





006 


90.0 ± 


0.6 


2 



Col. 1 : Cluster name. Col. 2: Cluster redshift. Cols. 3-6: fit parameters of the elliptical 2D p model: x c , y c is the central 
position; q is the projected axial ratio, and 9 is the orientation angle (north over east). In the last column, label 1 is 
for Chandra and 2 for XMM observations. 



tion of the sample clusters. The axial ratios are cal- 
culated with the inclination set to the median of the 
band of inclinations that are deemed suitable for a 
given cluster. 

7. Discussions 

In this work we have presented a simple but novel 
technique to extract the intrinsic shape and incli- 
nation class of galaxy clusters, using X-ray bright- 
ness maps alone. The availability of SZe data is then 
shown to improve the inclination estimate consider- 
ably. The main motivation behind this exercise is the 
improvement on mass estimates which is expected to 
refine the constraints that we can place on cosmolog- 
ical constants, from the analysis of cluster data. 

Our suggested formalism relies upon the inter- 
comparison of the density profiles that are depro- 
jected from the measured X-ray surface brightness 
maps, under 4 different deprojection scenarios that 
combine the extrema of the inclination scale and 
the geometry scale, when clusters are treated as fig- 
ures of revolution. It is found that each type of ob- 
served system leaves its unique signature in this 
inter-comparison amongst the deprojected X-ray lu- 
minosity density profiles. The knowledge of the clus- 



ter elongation along the LOS (from SZe data) is then 
used to supplement our shape and inclination deter- 
mination; in particular, the recovered constraints on 
the inclination are significantly improved, if we be- 
lieve the SZe measurements. In case of triaxial sys- 
tems, the intrinsic axial ratios can also be tracked, 
using the SZe data. 

7.1. Differe nces with Results of \Sereno et all 

I20M) 

Our shap e determination agre es with the shape ad- 
vanced bv lSereno et alj (|2006), for 13 of the 24 sam- 
ple clusters. This quantification i ncludes those cases 
that we identify as triaxial and ISereno et aL I (120061) 
find compatible with both the prolate and oblate ge- 
ometries. For another 10 cases, while we find the 
system to be definitely triaxial, (ISereno et al. 2006) 
report the system to be prolate . This might result 
from the fact that ISereno et al.l (120061) use the ob- 
served values of projected axial ratio and LOS extent, 
in equations relating these quantities and the intrin- 
sic axial ratio to inclination i, for the two geometries 
of prolateness and oblateness (see Equations 8-11 in 
Sereno et. al, 2006). 
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Table 2. Shapes & Inclinations of Cluster Sample 





Inferred Shape 


Approximate Incl 
(degrees) 


Shape in Sereno et al. (2006) 


Incl in Sereno et al. (2006) 
(degrees) 


MS 1137.5+6625 


prolate 


26 < / < 90 


prolate 


16+9 


MS 0451.6-0305 


triaxial-Type I 


40 < i < 90 


prolate/oblate 


56+13/58+43 


CI 0016+1609 


triaxial-Type I 


34 < i < 90 


prolate 


26±10 


RXJ1347.5-1145 


triaxial-Type III 


i < 47 


oblate/prolate 


84+117/35+12 


A 370 


triaxial-Type II 


i < 50 


prolate 


27±12 


MS 1358.4+6245 


prolate 


i < 41 


oblate 


52±8 


A 1995 


triaxial-Type III 


i < 36 


oblate/prolate 


64±35/46±27 


A611 


prolate 


29 < i < 90 


prolate 


35±32 


A 697 


triaxial-Type I 


42 < i < 90 


prolate 


26±10 


A 1835 


prolate 


35 < / < 90 


prolate 


29±7 


A 2261 


prolate 


i < 12 


prolate 


10±6 


A 773 


triaxial-Type II 


i< 36 


prolate 


14±5 


A 2163 


triaxial-Type I 


34 < / < 90 


prolate 


28±10 


A 1689 


prolate 


29 < / < 90 


oblate/prolate 


65±44/46±23 


A 665 


triaxial-Type I 


36 < i < 90 


oblate 


47±13 


A 2218 


triaxial-Type I 


34 < / < 90 


prolate 


22+11 


A 1413 


triaxial-Type II 


i < Al 


oblate/prolate 


73±33/37±15 


A 2142 


triaxial-Type II 


< 50 


oblate/prolate 


79±42/34±10 


A 478 


triaxial-Type II 


i < 48 


prolate 


23+12 


A 1651 


prolate 


i < 32 


prolate 


12+7 


A 401 


triaxial-Type II 


i < 40 


prolate 


25±6 


A 399 


oblate 


34 < i < 90 


oblate 


38+5 


A 2256 


triaxial-Type II 


i < 41 


prolate 


34±15 


A 1656 


prolate 


i < 29 


prolate 


33±30 



Col. 1: Cluster name, Col. 2: Intrinsic cluster shape recovered from our work, Col. 3: Recovered inclination class - 
in degrees, Col. 4: Shape reported by Sereno et. al (2006) & Col. 5: Inclination as reported by Sereno et. al (2006) 
(in degrees) - the clusters for which both prolate and oblate solutions were found admissible bv lSereno et al.l J2006I) 
are marked accordingly in Col 4 and the corresponding inferred inclinations are shown in Col 5. 



Table 3. Improved Inclinations (using SZe Data) of the Prolate and Oblate Sample Clusters 





Inferred Shape 


Incl from Deprojection 


from SZe data 


Refined Incl 






(degrees) 


(degrees) 


(degrees) 


MS 1137.5+6625 


prolate 


26 < / < 90 


1 1 < i < 63 


26 < i < 63 


MS 1358.4+6245 


prolate 


i < 41 


< i < 44 


41 < i < 44 


A611 


prolate 


29 < / < 90 


< i < 20.665 




A 2261 


prolate 


i < 12 


7 < i < 20 


7 < i < 12 


A 1689 


prolate 


29 < i < 90 


< i < 33 


29 < i < 33 


A 1651 


prolate 


i < 32 


8 < i < 32 


8 < i < 32 


A 1656 


prolate 


i < 29 


19 < / < 34 


19 < i < 29 


A 1835 


prolate 


35 < i < 90 


13 < i < 78 


35 < i < 78 


A 399 


oblate 


34 < i < 90 


34 < i < 37 


34 < i < 37 



Col. 1: Cluster name, Col. 2: Intrinsic cluster shape recovered from our work, Col. 3: Inclination recovered from 
the deprojection exercise, C ol. 4: Inclination as indicated by relating the projected POS axial ratio and LOS ex- 
tent dDe Filippis et alj |2005). Col. 5: The most improved bounds on inclination, obtained as overlap of constraints 
imposed by the deprojection exercise and the SZe data. 



They constrain their choice of geometry by iden- 
tifying realistic solutions for inclinations, such as 
— 1 < cos(/) < 1 . Thus, if for example, they identify 
an unrealistic solution for i for the oblate case, they 
declare the system prolate. But strictly speaking, all 
that they can conclude is that the cluster at hand is 
not oblate. In other words, it is possible that the clus- 
ter is triaxial, in such a way that given the high mea- 



surement errors, the equation for i in the prolate case 
also yields acceptable value s. The 10 sample clu sters 
that we spot as triaxial, but ISereno et alj d2006l) call 
prolate, are similarly 10 non-oblate clusters. 

The one s ystem that we identify as prolate and 
I Sereno et al.l d2006l) call oblate i.e. non-prolate, is 
MS 1358.4+6245. However, when we use the val- 
ues of 1.325 for the projected axial ratio (ignoring 



14 Chakrabarty, De Filippis & Russell: Cluster Geometry & Inclination 

Table 4. Axial Ratios and Improved Inclinations (using SZe Data) of the Triaxial Sample Clusters 





Shape from Deprojection 


Approx Incl from Deproj 
(degrees) 


Correct Incl 
(degrees) 


H 


92 


MS 0451.6-0305 


triaxial-Type I 


40 < 


< 90 


67 < / < 90 


0.92 


1.34 


CI 0016+1609 


triaxial-Type I 


34 < 


< 90 


73 < / < 90 


0.65 


1.24 


A 697 


triaxial-Type I 


42 < 


< 90 


77 < i < 90 


0.58 


1.33 


A 2163 


triaxial-Type I 


34 < 


<90 


72 < i < 90 


0.68 


1.24 


A 665* 


triaxial-Type I 


36 < 


<90 


70 < / < 90 


0.78 


1.27 


A 2218 


triaxial-Type I 


34 < 


<90 


78 < / < 90 


0.45 


1.20 


A 370 


triaxial-Type II 


50 < 


<90 


73 < i < 90 


0.84 


1.62 


A 773 


triaxial-Type II 


36 < 


<90 


78 < i < 90 


0.48 


1.28 


A 1413 


triaxial-Type II 


47 < 


<90 


66 < ( < 71 


0.96 


1.64 


A 2142 


triaxial-Type II 


50 < 


<90 


69 < / < 77 


0.96 


1.67 


A 478 


triaxial-Type II 


48 < 


<90 


75 < / < 90 


0.73 


1.52 


A 401 


triaxial-Type II 


40 < 


< 90 


72 < i < 90 


0.76 


1.34 


A 2256 


triaxial-Type II 


41 < 


<90 


67 < / < 90 


0.94 


1.36 


A 1995 


triaxial-Type III 


i< 36 


< i < 31 


3.94 


0.34 


RXJ1347.5-1145 


triaxial-Type III 


i < 47 


< i < 65 


5.14 


0.29 



Col. 1: Cluster name, Col. 2: Intrinsic cluster shape recovered from our work, Col. 3: Approximate inclination range, 
guessed from the deprojection exercise, Col. 4: Bounds on inclination, obtained from the constraints on the intrinsic 
axial ratios, as indicated by the type of triaxiality (Type I, II or III) that is assigned to the cluster from deprojection 
work, Col. 5: Ratio q t between the length of the principal axis (a) that corresponds to one of the photometric axes and 
that along the y-principal axis, (b, see Figure[3} Col. 6: Ratio (q 2 ) between a and and the other remaining principal 
axis (c). 

* For this cluster, the q L os value reported at the higher end of the range is used; the medial value of q L os constrains 
the cluster inclination to less than a degree. 



the +1% error) and 0.91 for e L0S (using the upper 
bound from the measured range of 0.72+0.19), in 
Equations 8 and 9, we do actually get a realistic in- 
clination of about 44°, under the assumption of pro- 
lateness. In other w ords, the methodology used by 
ISereno et al.l (2006), does not rule out prolateness for 
this cluster, in line with our inference. 

Thus, our methodology is relatively more pow- 
erful, since we have greater resolution ability than 
simply ruling out prolateness/oblateness. Moreover, 
our implementation of the SZe data allows for the re- 
covery of a much narrower range of inclinations, for 
the same measurement as used by previous workers. 

7.2. X-rays vs. Optical SB 

We would like to emphasise that it is very much pos- 
sible to utilise the optical surface brightness distri- 
bution of clusters, if available, to make an estimate 
of the morphological structure of the cluster with the 
aid of DOPING. This would allow the verification 
and quantification of possible discrepancies between 
results obtained using X-ray and optical data. This 
is particularly interesting, giv en the study presented 
by IGottlober & Yep es (2007) which suggests a rel- 
atively more spherical central gas component, com- 
pared to the dark matter one. 

To this aim, we are now applying the DOPING 
algorithm to both dark matter halos from the 
Millennium Simulations and to gas+dark matter SB 
distribution of the same halos (from the Millennium 



Gas Simulations). Results will be published in 
Chakrabarty et. al, (2008). This will answer the ques- 
tion of how the shape determination is affected by 
changes in the nature of the input. 

In this connection, it is further important to men- 
tion that our sample of real clusters show more ec- 
centric gas distribution in projection, than the simu- 
lat ed ones studied in the hy drodynamics simulations 
of IGottlober & Yepesl (|2007|) and, consequently, the 
intrinsic axial ratios that we recover, are not akin to 
#1 = 1, qi =2, as suggested in that work. 

7.3. Core Sizes 

We tested our predictions on a suite of model clusters 
and found that the recovery of shapes and inclina- 
tions is independent of the axial ratios that were as- 
signed to these test systems, as well as the core sizes 
that these test systems were described by. That the 
core size is not influential in the recovery of the in- 
trinsic shape of the cluster is not surprising since it is 
the outer part of the brightness profile that is conven- 
tionally considered when identifying shapes of clus- 
ters. Thus, whether it is Abell 2261 or MS 1358, our 
formalism is able to do justice to systems at both ends 
of core sizes in the sample used herein. This brings 
us to an important point in regard to the applicability 
of our technique - it can be implemented to estimate 
the 3-D geometry of clusters at varying redshifts! In 
fact, the sample that we used to illustrate the efficacy 
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of our formalism is rather eclectic in its redshift cov- 
erage, (z=0.023 to 0.784). 

Crucially, our work has clearly indicated the de- 
pendence of the recovered core size of the depro- 
jected luminosity density profile, on the deprojection 
scenario used. Thus, for Abell 665, the deprojected 
density corresponding to an assumption of oblate- 
ness and i = i m j n implies a core size that is less than 
half (* 0.06 arcmin) that obtained by deprojection 
performed under prolateness and i = i,„j„ and incli- 
nation combinations («0.15 arcmin). 

It is important to stress here that we do not need 
to delete the cores from our analysis due to any con- 
ceived inability of our code to deal with the mod- 
elling of the very central regions of the clusters. The 
code used in the work is in fact non-parametric and 
its functionality is not challenged by local changes 
in slopes of the observed X-ray brightness profiles 
or even by local enhancements of observational er- 
rors in the brightness distributions. The deprojected 
profile will of course bear larger error bars when the 
input brightness profile manifests the same, than oth- 
erwise. The fact that the inclusion of the core does 
not affect the shape determination is brought home 
by the tests that we have performed by varying the 
core sizes of toy clusters. 

7.4. Central Density 

One ramification of the observed nature of the 
luminosity-temperature relationship in clusters is the 
need to invoke some (unknown) pressure that is ef- 
fective in red ucing the density of the intracluster gas 
in the cores (IVoitet al.ll2002h . However, as we have 
seen in this work, the X-ray luminosity density at 
the centres of clusters depends critically on the de- 
projection scenario adopted in the model. Thus, in 
Abell 1689, we notice that the density recovered un- 
der an assumption of prolateness and i = i m in, at 
about 3 arcmin, is about 10 times less than that re- 
covered under the other three deprojection scenarios. 
Thus, our work indicates the important contribution 
that deprojection uncertainties can make towards the 
observed trends in the cluster self-similar relation- 
ships. 

7.5. Effect of Errors of Cluster Elongation 

The current state of affairs regarding the quality of 
SZe data is indeed unsatisfactory, albeit improve- 
ments are impending. Thus, when we use the SZe 
data in this work, we make the conscious decision to 
work with the value of the cluster elongation that is 
reported at the centre of the error band. Incorporation 
of the measurement errors would have greatly re- 
duced the quality of the constraints on the inclina- 
tions. 



7.6. The Case of Abell 1995 and RXJ1347 

The clusters Abell 1995 and RXJ1347 have been 
found to be triaxial, of Type III, which implies that 
these clusters are nearly prolate except that c < a; 
also, we expect b - a <K a. The density profiles re- 
covered under the assumption of prolateness are flat- 
ter for these clusters than the profiles deprojected un- 
der oblateness. This comparison is allowed if the el- 
lipsoidal radius in the prolate cases, in general, ex- 
ceed that in the oblate cases. This is in turn ensured 
if a is the photometric semi-axis and if the inclina- 
tion is small. We find that the inclination is indeed 
constrained to very small angles, given the reported 
values of q p and qios (rather the medial value of 
these quantities within the reported measurement er- 
ror bands). Thus, these two clusters are deciphered to 
be nearly face-on systems. 

7.7. The Case of Abell 370 

The cluster Abell 370 was reporte d to bear a "pro- 
nounc ed triaxial morphology" by [De Filippis et al.l 
(2005), o n the basis o f the X-ray and SZe data. 
However, ISereno et al.l (120061) find the prolate solu- 
tion to be consistent with the observed X-ray and SZe 
data for this cluster. In contrast to this, we actually 
find this cluster to be triaxial, with an inclination in 
excess of 50°. 

7.8. The Case of CL 0016+609 

iHughes & Birkinshawl dl998l) suggested a "rea- 
sonable triaxial" morphology for the cluster 
CL 0016+609, on the basis of the discrepancy 
between the Hubble constant values that were 
estimated from the oblate and prolate models. This 
suggestion is in line with our conclusion of a triaxial 
shape for this cluster. 

7.9. Future Work 

The identification of the intrinsic shape and inclina- 
tion of a cluster is a major step in the characterisation 
of clusters, both in terms of cluster masses and quan- 
tification of the contribution of deprojection uncer- 
tainties to the observed scatter in the self-similar rela- 
tionships of clusters. These applications are planned 
in future contributions. 

In fact, it is envisaged that the improved under- 
standing of cluster spatial configurations leads to bet- 
ter constraints on cluster masses. If this is supple- 
mented with information of mass distributions from 
dynamical considerations (using an algorithm such 
as CHASSIS - Chakrabarty & Saha, 2001), we could 
potentially place bounds on the distribution of gas 
that is at hydrostatic equilibrium. Furthermore, using 
estimation of the mass within a given projected ra- 
dius from lensing measurements, would constrain the 
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dynamical mass distribution even better, leading to 
added improvements in cluster characterisation. An 
exercise in the determination of improvement in the 
dynamical mass of a lensing galaxy, using such lens- 
ing constraints, is currently underway (Chakrabarty, 
et. al, 2008). 
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Appendix A: 

Given that in our coordinate system, X — x, we say: 



(y,z) = Ri(Y,Z) 



(A.l) 



where R, is the rotational matrix corresponding to a 
rotation through the angle i. Let us consider a triaxial 
system, with intrinsic axial ratios q\ and <?2, i.e. ratios 
of the maximum extent along the X-axis to that along 
the F-axis is q\ and the ratio of the principle axis 
along the X-axis to that along the Z-axis is qi- 

In this triaxial configuration, the square of the el- 
lipsoidal radius is 



£ = x 2 + y 2 [q\ cos 2 (i) + q\ sin 2 (7)] 



(A.2) 



+Z 2 [q\ sin 2 (0 + q\ cos 2 (/)] + yz sin(2;')(<7 2 - q\) 

This implies that the plane-of-the-sky projection 
I(x,y) of any intrinsic quantity G(£), can be repre- 
sented as: 



.y) = 

Jo 

-f 

J, I 



CK&dz 



(A3) 



G(0 



d£ 



where 



2 2 

2 , 2^1^2 

v = x +y 



(A.4) 



The notation here is chosen to concur with that in 
iFabricant etall (11984 . Since the RHS of this equa- 
tion is a function of rj only and the projected quantity 
I(x,y) is a constant along the isophotes, the equa- 
tion relating rj to x and y must give the form of the 
isophotes. Then the ratio q p , of the extent along the 
x-axis to that along the >'-axis is related to ratio of the 
semi-axes along X, K&Z-axes respectively as: 



4 



q^ cos 2 (/) + q 2 sin (/) 



(A.5) 



A similar exercise involving integration over y yields 
the ratio e^os between the extent along the z-axis and 
the jc-axis as: 



l l 



,2 

'LOS 



cos 2 (/) sin (i) 

5 ' 5 



(A.6) 



The cluster elongation index qws that is used 
elsewhere in the pape r (notation borrowed from 
iDe Filippis et al.l (120051) ) differs from e L0S in that 
qws is the ratio along the LOS to the photometric 
major axis. 

Thus, for oblate systems, when q\ — 1, qi — q > 
1, q p > 1 and q LO s = e ws , 



<r 



q 2 cos 2 (/) + sin (?) 



(A.7) 
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Using this in Equation IA.6I we get that for oblate 
systems, 



1 - q\ cos 2 (/) 
sin 2 (/) 



qios = (a.8) 



For prolate systems, q\ = 1, qi < 1 ==> qi = 
l/<7, q P < 1 and q LO s = e LO s/e p . So, 

ql = cos 2 (/) + q 2 sin 2 (/) (A.9) 



Using this in Equation IA.6I we get that for prolate 

systems, 

q 2 n - COS 2 (i) 

ilos = ^V-vT (A ' 10) 

cf p sin (0 

Since q p and qios are observables, q\ and qi can 
be calculated for the three different types of triaxial 
systems, as long as we keep in mind that: 

- for Triaxial-Type I: q\ < \,qi > 1, q p < 1 and 

qws = ews lip- 

- for Triaxial-Type EE: q\ < 1, qi > 1, q p > 1 and 

<7los = e ws ■ 

- for Triaxial-Type HI: #i > 1, qi < 1, q p > 1 and 



We use Equations IA.5I and IA.6I to determine q\ 
and q2 in a triaxial configuration. To do this, the val- 
ues of q p and qios will of course need to be sup- 
plied from observations. The value of qios that we 
consider is the centroid of the reported error band in 
iDe Filippis etail (120051) . 

However, for the different types of triaxial sys- 
tems that we consider, there are two other constraints 
involving q\ and q2 (in the form of inequalities); 
such inequalities are typical of the triaxial type of 
the cluster at hand and are listed above for Types I, II 
and III. Implementing these would imply two distinct 
ranges of inclinations for each cluster. The overlap 
of these provides the final constraints on the incli- 
nations. Once this final range of inclinations is ob- 
tained, Equations IA.5I and IA.6I are solved to give q\ 
and qi, at the median of the recovered inclination 
range. 



